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ABSTRACT:  This  project  serves  as  a  demonstration  of  the  three-dimensional  (3D)  numerical  modeling 
capability  of  large  reservoirs  up  to  and  including  the  complex  geometries  of  the  powerhouse,  spillway, 
and  lock  guard  wall  structures.  The  3D  Navier-Stokes  module  of  the  Adaptive  Hydraulic  (ADH)  computer 
code  is  used  to  compute  the  forebay  flow  at  the  Ice  Harbor  Dam.  A  numerical  flow  model  of  the  Ice 
Harbor  Dam  forebay,  which  is  located  on  the  Snake  River,  was  constructed.  The  model  results  are 
compared  to  field  data  and  solutions  obtained  by  an  independent  laboratory’s  in-house  code. 
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Preface 


The  model  investigation  reported  herein  was  authorized  and  funded  as  a 
demonstration  project  for  the  Navigation  Systems  Research  Program.  Additional 
funding  for  model  validation  was  provided  by  the  U.S.  Army  Engineer  Walla 
Walla  District  as  a  part  of  the  Lower  Monumental  Dam  forebay  computational 
fluid  dynamics  (CFD)  study  (in  progress). 

This  study  was  conducted  at  the  Coastal  and  Hydraulics  Laboratory  (CHL), 
U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  during  the 
period  November  2003  to  March  2005  by  Drs.  Richard  L.  Stockstill,  John  E. 
Hite,  Jr.,  R.  Charlie  Berger,  Mr.  Alex  R.  Carrillo,  and  Ms.  Jane  M.  Vaughan. 

Dr.  Larry  Weber,  Director  IIHR-Hydroscience  and  Engineering,  University  of 
Iowa,  and  Dr.  Songheng  Li,  Assistant  Research  Engineer,  IIHR-Hydroscience 
and  Engineering,  University  of  Iowa,  served  as  independent  technical  reviewers. 
Dr.  Sandra  K.  Knight  was  the  Technical  Director  for  navigation  and  Mr.  James 
E.  Clausner  was  the  manager  of  the  Navigation  Systems  Research  Program. 

The  work  was  performed  under  the  general  supervision  of  Mr.  Thomas  W. 
Richardson,  Director,  CHL;  Dr.  William  D.  Martin,  Deputy  Director,  CHL;  Dr. 
Rose  M.  Kress,  Chief  of  the  Navigation  Division,  CHL;  Mr.  Bruce  A.  Ebersole, 
Chief  of  the  Flood  and  Storm  Protection  Division,  CHL;  Mr.  Donald  C.  Wilson, 
Chief  of  the  Navigation  Branch,  CHL;  and  Dr.  Robert  T.  McAdory,  Chief  of  the 
Estuarine  Engineering  Branch,  CHL.  Drs.  Stockstill  and  Hite  and  Ms.  Vaughan 
wrote  this  report. 

At  the  time  of  publication  of  this  report  Dr.  James  R.  Houston  was  Director 
of  ERDC,  and  COL  James  R.  Rowan,  EN,  was  Commander  and  Executive 
Director. 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

miles  (U.S.  statute) 

1.609344 

kilometers 

square  feet 

0.09290304 

square  meters 

cubic  feet 

0.028316847 

cubic  meters 
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Introduction 


This  investigation  serves  as  a  demonstration  of  the  three-dimensional  (3D) 
numerical  modeling  capability  of  large  reservoirs  up  to  and  including  the  com¬ 
plex  geometries  of  the  powerhouse,  spillway,  and  lock  guard  wall  structures.  The 
3D  Navier-Stokes  module  of  the  ADaptive  Hydraulic  (ADH)  computer  code 
(Stockstill  and  Berger  2001,  Berger  and  Stockstill  1999)  is  used  to  compute  the 
forebay  flow.  ADH  is  a  suite  of  finite  element  models  of  which  the  Navier- 
Stokes  module  solves  the  Reynolds-Averaged  Navier-Stokes  equations  using  the 
method  of  finite  elements.  The  solution  provides  the  three  components  of  veloc¬ 
ity  and  the  pressure  at  each  node  in  the  computational  mesh.  Secondary  variables 
such  as  velocity  gradients  can  also  be  computed  for  each  node. 

Ice  Harbor  Dam  forebay  (Washington)  was  selected  as  the  project  on  which 
the  modeling  capabilities  were  to  be  demonstrated  because  velocity  data 
throughout  the  forebay  were  available  for  model  validation.  The  Ice  Harbor  Dam 
is  located  on  the  Snake  River  about  9.7  miles  upstream  of  its  confluence  with  the 
Columbia  River  as  shown  on  the  vicinity  map  on  Figure  1 .  Flow  conditions  at  the 
project  while  spilling  are  shown  in  Figures  2  and  3. 


Project  Features 

The  Ice  Harbor  project  was  authorized  by  Section  2  of  the  River  and  Harbor 
Act  of  1945  and  was  approved  March  2,  1945.  The  lake  behind  Ice  Harbor  Dam 
has  been  designated  Lake  Sacajawea.  Ice  Harbor  Lock  and  Dam-Lake  Sacajawea 
is  Unit  1  of  4  included  in  the  Lower  Snake  River  Project,  Washington,  and  Idaho. 
The  dam  is  2,822  ft1  long,  with  an  effective  height  of  100  ft.  It  is  of  the  concrete 
gravity  type,  with  an  earthfill  embankment  section  at  the  north  abutment.  The 
dam  has  a  ten-bay  spillway  that  is  590  ft  long,  and  ten  50-ft  tainter  gates.  At 
el  440, 2  Lake  Sacajawea  extends  northeast  32  miles  upstream  to  Lower 
Monumental  Dam.  The  Ice  Harbor  powerhouse  contained  three  90,000-kW  units 
initially.  It  now  has  three  additional  1 1 1,000-kW  units  and  is  capable  of  gener¬ 
ating  603,000  kW  of  power.  The  lock  has  clear  chamber  dimensions  of  86  ft  by 
675  ft,  and  a  103-ft  lift.  The  structure  has  two  fish  ladders  for  passing  migratory 
fish.  Construction  of  the  Ice  Harbor  Project  began  in  December  1955,  and  the 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on  page  vi. 

2  All  elevations  (el)  cited  herein  are  in  feet  referenced  to  the  National  Geodetic  Vertical  Datum 
(NGVD). 
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project  began  operation  in  December  1961.  Powerhouse  units  4  though  6  were 
installed,  with  all  units  producing  power  by  January  1976. 


Purpose  and  Scope 

The  purpose  of  this  study  was  to  evaluate  ADH’s  ability  to  simulate  forebay 
flow.  A  numerical  flow  model  of  the  forebay  approaching  the  Ice  Harbor  Dam 
was  constructed  and  two  project  flow  conditions  were  simulated.  In  particular, 
velocity  comparisons  were  made  with  field  data  supplied  by  the  U.S.  Army 
Engineer  District,  Walla  Walla,  office.  The  study  also  tested  the  model’s  sensi¬ 
tivity  to  selection  of  eddy  viscosity  and  grid  resolution.  Secondary  flow  variables 
of  acceleration  and  velocity  gradients  were  then  compared  with  those  computed 
by  another  flow  solver. 

This  report  is  a  summary  of  the  Ice  Harbor  Dam  forebay  computational  flow 
model  and  serves  as  a  demonstration  of  the  capability  to  simulate  the  flow 
approaching,  engaging,  and  passing  the  project. 
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2  Model  Description 


The  computational  model  consists  of  the  governing  equations,  the  discretiza¬ 
tion  scheme  used  to  numerically  solve  the  equations,  the  computational  mesh  on 
which  the  domain  is  discretized,  and  the  boundary  and  initial  conditions  needed 
to  close  the  system  of  equations. 


Governing  Equations 

The  Reynolds-Averaged  Navier-Stokes  (RANS)  equations  are  employed  to 
model  the  flow  field  approaching,  interacting,  and  passing  hydraulic  structures. 
The  RANS  equations  are  three-dimensional  (3D)  with  four  degrees  of  freedom: 
the  pressure  and  the  three  components  of  fluid  velocity.  The  RANS  equations 
make  no  assumptions  as  to  pressure  distributions  and,  because  many  hydraulic 
flow  models  assume  the  flow  is  hydrostatic,  RANS  models  are  referred  to  as  non¬ 
hydrostatic  models. 

A  finite  element  discretization  is  employed  to  solve  the  RANS  equations. 

The  particular  code  is  the  Navier-Stokes  module  of  the  ADaptive  Hydraulics 
(ADH)  flow  model  developed  in  the  Coastal  and  Hydraulics  Laboratory  (CHL) 
of  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC).  Tetrahe¬ 
dral  shaped  elements  are  used  in  the  finite  element  formulation.  This  unstructured 
approach  is  quite  flexible  in  representing  the  flow  field’s  geometry. 

If  the  instantaneous  velocity  is  split  into  a  mean,  V,  and  a  fluctuating  compo¬ 
nent,  u ,  the  incompressible  Navier-Stokes  equations  can  be  time  averaged  to  pro¬ 
duce  the  RANS  equations.  The  RANS  equations  are  written  in  terms  of  the  mean 
velocity,  V(x,t),  and  pressure,  p(x,t),  so  as  to  reduce  the  modeling  of  turbulence  to 
a  set  of  equations  that  incorporate  terms  to  model  the  effects  of  turbulence  on  the 
mean  flow.  The  mean  and  fluctuating  velocity  and  pressure  are  then  substituted 
into  the  mass  and  momentum  equations,  and  after  time  averaging  over  intervals 
that  are  long  compared  with  the  periods  of  fluctuations,  one  obtains  the  RANS 
equations: 


V  -V  =  0 


(i) 


and 
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where 


x  =  spatial  coordinate 
t  =  time 

p  =  fluid  density 
a  =  -/;/+t 
r  2pD 

D=  ^(Vu  +  VuT) 

I  =  identity  matrix 
//  =  fluid  viscosity 

In  a  RANS  approach,  the  term  V  •  ( puu )  is  used  to  represent  the  effect  of  tur¬ 
bulence  on  the  mean  flow.  Following  the  suggestion  of  Boussinesq,  an  eddy  vis¬ 
cosity  is  added  to  the  molecular  viscosity  in  the  momentum  equations  to  account 
for  the  effects  of  turbulence.  An  algebraic  eddy  viscosity  model  (also  known  as  a 
zero  equation  model)  was  used  for  the  ADH  demonstration  application  to  the  Ice 
Harbor  Dam  forebay.  ADH  also  has  a  k-s  turbulence  closure  model,  but  this  more 
sophisticated  model  was  not  used  for  this  demonstration.  The  k-s  model,  in  par¬ 
ticular,  calculates  the  turbulent  kinetic  energy,  k ,  and  its  dissipation  rate,  s , 
throughout  the  domain.  These  two  variables  are  used  to  define  an  eddy  viscosity, 
pt.  The  Reynolds  stress,  r#,  can  be  expressed  using  the  kinematic  form  of  the  eddy 
viscosity. 


Physical  Description 

Before  the  equations  of  motion  can  be  applied,  the  domain  must  be  discre¬ 
tized.  This  process  includes  the  construction  of  a  3D  computer-aided  drafting 
(CAD)  representation  of  the  flow  boundaries  including  the  geometric  features  of 
the  hydraulic  structure  and  the  bathymetry  of  the  approaching  river.  The  3D  sur¬ 
face  model  of  the  hydraulic  structure  and  the  forebay  bathymetry  CAD  model 
was  constructed  using  AutoCAD1  (http://usa.autodesk.com/).  The  coordinates 
were  in  a  state  plane  coordinate  system  where  the  x-direction  is  easterly  and  the 
y-direction  is  northern.  The  entire  model  was  translated  2,410,000  ft  toward  the 
west  (i.e.  Ax  =  -2,410,000  ft)  and  330,000  ft  toward  the  south  (i.e.  Ay  =  -330,000 
ft).  This  provides  more  significant  digits  and  can  produce  more  accurate  solu¬ 
tions.  The  CAD  model  was  then  used  as  input  for  the  grid  generator.  The  3D  vol¬ 
ume  mesh  is  of  the  unstructured  type  using  tetrahedral  elements.  Figures  4-10  are 
pictures  of  particular  portions  of  the  Ice  Harbor  CAD  model.  Approximately  one 
mile  of  the  Snake  River  upstream  of  the  dam  was  reproduced.  The  river  is  a  little 


1  Statement  of  product  names  does  not  infer  product  endorsement. 
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over  2,000  ft  wide,  immediately  upstream  of  the  lock  guard  wall.  The  model 
reproduced  the  10  spillway  bays,  complete  with  piers  and  tainter  gates  in  the 
appropriately  raised  position,  and  the  6  powerhouse  units  each  having  three 
intakes  to  a  particular  turbine.  The  trash  racks  were  not  included  in  the  model, 
but  the  features  such  as  the  slots  and  other  important  offsets  were.  The  navigation 
lock’s  floating  guard  wall  protrudes  into  the  flow  field  and  was  included  in  the 
model. 

The  spillway  portion  of  the  model  included  the  face  and  crest,  up  to  the  point 
where  the  flow  exited  under  the  opened  tainter  gates.  The  piers,  walls,  and  soffits 
of  the  powerhouse  intakes  were  included.  Forty-two  feet  of  the  conduit,  from  the 
powerhouse  face  to  the  end  of  the  soffit’s  curved  surface,  was  reproduced. 


Computational  Mesh 

A  computational  mesh  consisting  of  609,721  tetrahedral  elements  and 
153,602  nodes  was  constructed  to  fill  the  volume  created  by  the  CAD  model.  The 
commercial  mesh  generator  used  for  the  Ice  Harbor  project  was  ANSYS  ICEM 
CFD  (http://www-berkeley.ansys.com/).  Element  sizes  varied;  side  lengths  ranged 
from  nearly  30  ft  near  the  center  of  the  river  about  2,500  ft  upstream  of  the  dam 
to  a  few  inches  at  the  intake  to  the  powerhouse  and  under  the  tainter  gates  on  the 
spillway.  The  surface  meshes  shown  in  Figures  11-15  partially  illustrate  the 
resolution  of  the  computational  mesh.  The  surface  mesh  is  composed  of  the  indi¬ 
vidual  faces  of  the  tetrahedral  elements  that  form  the  boundary.  Boundary  condi¬ 
tions  of  these  faces  and  their  nodes  are  needed  to  determine  a  particular  solution 
to  the  governing  partial  differential  (RANS)  equations.  Natural  boundary  condi¬ 
tions  involving  fluxes  are  assigned  to  the  surface  faces  and  Dirichlet  boundary 
conditions  are  assigned  to  the  surface  nodes.  The  surface  mesh  and  portions  of 
the  volume  mesh  in  regions  near  the  dam  are  shown  in  Figures  16  and  17. 


Computational  Issues 

The  ADH  model  solves  the  time-dependent  RANS  equations.  Steady  state 
solutions  are  obtained  by  starting  from  some  assumed  initial  conditions  (usually  a 
quiescent  pool)  and  time  stepping  until  the  solution  reaches  steady  state  (if  a 
steady-state  solution  exists).  Often,  flow  fields  exhibit  shedding  eddies  and  such, 
and  there  is  no  steady-state  solution.  That  was  not  the  case  with  the  Ice  Harbor 
model.  The  Ice  Harbor  model  was  initiated  with  a  quiescent  pool,  and  given  the 
flux  boundary  conditions,  the  model  was  run  for  a  simulated  time  totaling  8 
hours.  At  this  point,  the  differences  in  velocities  and  pressures,  between  con¬ 
secutive  solutions,  were  small  enough  that  the  model  was  considered  converged 
to  steady  state. 


Computer 

The  computer  runs  for  this  project  were  primarily  carried  out  using  an  SGI 
Origin  3000  (03K)  Complex  housed  at  the  U.S.  Army  Engineer  Research  and 
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Development  Center,  Major  Shared  Resource  Center,  Vicksburg  site.  The  03K 
Complex  is  a  shared-memory  multiprocessor  system.  The  ADH  model  is  more 
universal  in  that  it  uses  distributed  memory  and  does  not  require  a  shared- 
memory  system,  but  the  03 K  was  readily  available  and  so  it  was  used  for  this 
project.  The  03K  Complex  consists  of  512  400-MHz  and  1,024  700-MHz  proc¬ 
essors  and  a  computational  capacity  of  1,833  GFLOPS.  Other  specifications  for 
this  computer  can  be  found  at  http:Zwww.wes. hpc.mil/index.htm. 


Partitioning 

Large  models  often  use  high-performance  computers  for  flow  solution,  and 
this  generally  involves  parallelization  of  the  problem.  The  method  of  dividing  the 
original  grid  among  many  processors  is  not  straightforward,  but  clever  division 
can  reduce  both  computation  time  and  communication  time  between  the  proces¬ 
sors.  The  basic  idea  is  to  divide  the  problem  in  a  manner  that  balances  the  loads 
applied  to  each  processor.  This  way,  all  processors  complete  a  particular  task  at 
about  the  same  time,  which  avoids  processors  sitting  idle  waiting  for  other  proc¬ 
essors  to  complete  their  assigned  tasks.  ParMETIS  ( http://www.users.cs.umn . 
edu/~karypis/metis/parmetis/index.html )  is  used  by  ADH  to  partition  the  mesh  in 
an  efficient  manner.  This  partitioning  is  dynamic  in  that  it  adjusts  as  needed  dur¬ 
ing  a  simulation  run.  Computer  runs  on  the  Ice  Harbor  forebay  model  generally 
employed  between  64  and  96  processors,  depending  on  the  current  load  on  the 
03K. 

Output 

From  a  practical  sense,  the  primary  variables  of  pressure  (scalar)  and  velocity 
(vector)  were  written  as  model  output  for  post-processing.  Other  computed 
quantities  such  as  the  error  indicator  used  for  the  adaption  trigger  (scalar)  and  the 
velocity  gradients  in  the  strain-rate  tensor,  important  information  for  the  Numeri¬ 
cal  Fish  Surrogate  (NFS),  were  also  written  to  files.  The  adaption  error  indicators 
are  useful  only  for  the  setup  of  mesh  adaption  and  the  strain  tensor  components 
were  stored  for  future  integration  with  the  NFS.  Sometimes,  it  is  informative  to 
plot  Froude  or  Reynolds  numbers,  and  these  values  can  be  computed  during  post¬ 
processing  procedure.  Solution  figures  for  this  report  were  generated  using 
FieldView  (http://www.ilight.com/)  and  Tecplot  (http://www.tecplot.com). 
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3  Validation 


Flow  Conditions 

As  a  demonstration  of  the  model’s  capability,  one  particular  flow  condition 
was  simulated.  The  particular  flow  condition  chosen  to  simulate  is  one  that  has 
been  successfully  modeled  by  the  IIHR-Hydroscience  and  Engineering  (IIHR). 
This  flow  condition  was  simulated  so  that  results  with  the  IIHR  model  results 
could  be  used  to  evaluate  the  ADH  simulations  as  to  their  reasonableness.  The 
pool  elevation,  tainter  gate  openings,  and  powerhouse  loadings  drive  the  dis¬ 
charges  listed  in  Table  1. 


Table  1 

Flow  Conditions  Simulated  for  Validation  with  Field  Data 

Total  River  Discharge  105,500  cfs,  Pool  el  439.5 

Powerhouse  Unit  (Left  to  Right 
Looking  D/S) 

Discharge, 

cfs 

Spillway  Bay  (Left  to  Right 
Looking  D/S) 

Discharge, 

cfs 

1 

14,900 

1 

1,700 

2 

15,400 

2 

3,500 

3 

0 

3 

3,500 

4 

16,500 

4 

5,200 

5 

0 

5 

5,200 

6 

18,700 

6 

5,200 

7 

5,200 

8 

4,000 

9 

3,900 

10 

2,600 

Boundary  Conditions 

The  boundary  conditions  employed  to  model  this  particular  flow  condition 
consist  of  setting  flow  rates  at  each  of  the  powerhouse  intakes,  each  spillway 
gate,  at  the  river  cross  section  making  up  the  upstream  limit  of  the  model,  at  the 
water  surface,  and  at  all  solid  boundaries  comprising  the  computational  domain. 
Each  spillway  bay  also  has  a  geometric  configuration  associated  with  the  opening 
of  the  bay’s  tainter  gate.  Boundary  conditions  at  the  solid  boundaries  prohibit  the 
flux  of  mass  and  momentum  through  the  boundary.  This  no-flux  boundary 
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condition  is  applied  to  the  solid  boundaries  and  a  drag  coefficient  is  assigned  to 
the  faces  of  each  element  forming  the  solid  boundaries.  The  water  surface  is 
modeled  as  a  rigid  lid,  whereby  it  is  fixed  in  space.  This  requires  that  the  pressure 
at  the  water  surface  be  calculated.  Modeling  the  water  surface  as  a  rigid  lid  is 
reasonable  for  large,  deep  bodies  of  water  where  the  Froude  number  is  small, 
such  as  the  case  of  the  flow  approaching  Ice  Harbor  Dam.  A  stress  free  condition 
is  applied  for  the  tangential  component  of  shear  at  the  water  surface. 

The  rigid  lid  assumption  seems  reasonable  for  the  problem  at  hand.  Here,  the 
water  surface  is  flat  in  part,  because  the  outflow  boundaries  are  all  submerged.  If 
the  simulation  were  to  include  an  outflow  at  or  near  the  water  surface  or,  if  for 
any  other  reason,  the  water  surface  is  expected  to  be  contorted,  then  the  Free  Sur¬ 
face  routines  in  ADH  could  be  employed.  ADH  uses  Space-Time  finite  elements, 
so  the  moving  mesh  routines  used  to  model  free-surface  problems  are  similar  to 
an  Arbitrary  Lagrangian-Eulerian  (ALE)  scheme  where  the  computational  mesh 
moves  such  that  free  surface  nodes  maintain  a  specified  pressure  (typically  zero). 
The  mesh  is  adjusted  each  time  step  in  response  to  the  surface  displacement.  This 
inhibits  element  shearing. 


Inflow  boundary  conditions 

Flux  boundaries  are  prescribed  as  shown  in  Table  1.  The  inflow  boundary  is 
the  sum  of  the  discharges  through  the  powerhouse  and  the  spillway.  In  the 
absence  of  specific  data  otherwise,  a  uniform  velocity  was  assumed  across  the 
inflow  boundary  as  the  total  discharge  divided  by  the  inflow  cross-sectional  area. 
This  uniform  distribution  was  used,  although  if  a  particular  distribution  is  known 
to  exist,  this  inflow  boundary  condition  could  be  setup  as  well. 


Outflow  boundary  conditions 

Outflow  discharges  from  the  hydraulic  features  listed  in  Table  1  were  speci¬ 
fied  as  natural  boundary  conditions  to  each  element  face  making  up  each  outflow 
boundary.  The  average  velocity  at  each  feature  was  applied  to  each  outflow  face. 

Specifying  both  inflow  and  outflow  mass  fluxes  in  a  fixed  domain  (rigid  lid), 
can  result  in  over  specification  of  boundary  conditions  because  of  roundoff 
errors).  To  avoid  over  specification  of  flux  boundary  conditions,  the  inflow 
boundary  condition  was  changed,  soon  after  model  startup,  from  a  specified  flow 
rate  to  a  simple  flux  boundary. 


Results 

General  results  of  the  flow  simulation  are  given  in  a  set  of  figures  of  velocity 
contours  and  vectors  at  horizontal  planes  sliced  through  the  flow  field  at  various 
elevations.  The  velocity  contours  for  each  figure  (Figures  18-37)  range  from  a 
minimum  value  of  0.0  fps  to  a  maximum  value  of  2.0  fps.  The  values  of  the  con¬ 
tour  lines  represent  the  numbers  provided  on  the  contour  scale.  Lower  velocity 
magnitudes  are  present  in  the  forebay,  especially  near  the  shallow  shores.  The 
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velocity  approaches  30  to  40  fps  as  the  flow  exits  under  the  spillway  tainter 
gates.  However,  for  visualization  purposes,  the  range  of  0.0  to  2.0  fps  seems  rea¬ 
sonable.  Figure  38  illustrates  the  flow  conditions  along  a  vertical  plane  placed 
normal  (perpendicular)  to  the  axis  of  the  dam  at  the  center  intake  of  powerhouse 
unit  6.  Figure  39  shows  the  velocity  distribution  on  a  vertical  plane  sliced  normal 
to  the  axis  of  the  dam  through  the  center  of  spillway  bay  6. 

Although  uniform  velocities  were  set  as  inflow  boundary  conditions,  the 
resulting  velocities  just  downstream  of  the  inflow  boundary  were  changed  dra¬ 
matically.  The  velocity  in  the  shallows  was  rapidly  reduced  from  those  imposed 
on  at  the  inflow  boundary.  The  radius  of  influence  is  quite  small  compared  to  the 
length  of  river  reproduced.  Modeling  of  other  projects  in  the  future  are  planned 
that  will  try  introducing  an  artificial  “box”  upstream  of  the  actual  topography. 
The  bed  elevations  from  the  flat  bottom  of  the  inflow  box  will  be  gradually  tran¬ 
sitioned  to  those  of  the  actual  project  bathymetry.  This  is  the  approach  used  in 
most  large  physical  models  and  should  prove  useful  in  a  computational  model  as 
well.  Then  a  uniform  velocity  of  total  discharge  divided  by  the  cross  sectional 
shape  of  the  inflow  “box”  will  be  more  appropriate  and  will  not  require  any  field 
data  by  which  the  inflow  boundary  conditions  are  to  adjusted. 

Figures  40  and  41  show  streamlines  of  flow  approaching  the  structure.  The 
powerhouse  units  whose  outflow  is  black  are  those  units  not  operating  (see 
Table  1). 


Comparison  of  Results 

A  comparison  of  the  CFD  results  was  made  with  prototype  data  at  selected 
locations  in  the  forebay.  The  electronic  files  of  the  prototype  data  at  depths  of 
approximately  10,  20,  40,  60  and  80  ft  below  the  water  surface  were  obtained  in 
the  form  of  spreadsheets  from  the  Walla  Walla  District.  The  water-surface  eleva¬ 
tion  was  el  439.5  during  collection  of  this  data.  The  discharge  through  the  spill¬ 
way  was  40,000  cfs  and  the  discharge  through  the  powerhouse  was  65,500  cfs. 

The  comparison  of  the  CFD  data  and  the  prototype  data  are  provided  in  Fig¬ 
ures  42-46.  The  red  vectors  are  the  prototype  data  and  the  black  vectors  are  the 
CFD  results.  Magnifying  the  screen  image  makes  the  comparison  easier  to  view. 
Figure  47  presents  a  typical  plot  of  the  prototype  data  without  the  CFD  results.  In 
general,  the  comparisons  are  favorable  and  show  the  CFD  model  does  a  good  job 
of  computing  the  velocity  vectors.  The  CFD  results  are  more  uniform  than  the 
field  data.  The  prototype  measurements  are  essentially  an  instantaneous  velocity 
and  are  not  time-averaged  over  the  turbulence  time  scales  as  those  of  the  CFD 
values.  The  direction  and  magnitude  of  the  prototype  data  such  as  that  shown  in 
Figure  47  appear  more  random  and  some  vectors  are  definitely  indicating  the 
wrong  flow  direction.  This  is  due  to  bad  samples  or  noise  in  the  data.  The  adja¬ 
cent  vectors  that  have  noticeably  different  directions  and  magnitudes  demonstrate 
the  unsteady  nature  of  the  flow  and  the  data  collection  process.  The  CFD  results, 
in  general,  under  predict  the  velocity  magnitude.  However,  the  comparison 
shows  that  the  CFD  solution  is  acceptable  for  making  engineering  decisions. 
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Grid  Resolution 


A  test  to  determine  the  effects  of  grid  refinement  was  conducted.  As  previ¬ 
ously  mentioned,  the  original  mesh  (hereto  referred  to  as  the  unrefined  grid)  was 
composed  of  153,602  nodes.  The  adaption  capabilities  of  the  ADH  model  were 
set  such  that  each  tetrahedron  in  the  volume  was  refined  one  level.  This  increase 
in  resolution  produced  a  mesh  containing  548,464  nodes.  Figures  48-50  show  a 
comparison  of  the  CFD  results  with  and  without  additional  grid  adaption  for 
depths  of  10,  40,  and  80  ft.  The  additional  grid  adaption  is  noted  as  the  refined 
CFD  results.  Figures  51  and  52  show  comparisons  of  velocity  magnitude  at  pow¬ 
erhouse  unit  6  and  spillway  bay  5,  respectively.  There  is  little  difference  in  the 
velocity  magnitudes,  though  a  slight  difference  is  evident  near  the  structure. 
Comparison  of  solutions  computed  on  the  original  mesh  with  those  computed  on 
the  refined  mesh  demonstrates  that  the  original  mesh  resolution  (153,602  nodes 
and  609,721  tetrahedra)  is  acceptable  for  computing  this  type  of  flow  field. 
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4  Model  Testing 


Flow  Conditions 

The  IIHR-Hydroscience  and  Engineering  has  been  successfully  modeling  the 
flow  in  forebays  of  projects  on  the  Lower  Snake  River  (i.e.,  Blank  and  Weber 
2000).  The  IIHR  model,  U2RANS  (Unsteady  and  Unstructured  Reynolds- 
Averaged  Navier-Stokes  solver),  has  been  validated  with  field  data  and  results 
from  U2RANS  have  been  compared  with  results  from  other  commercial  CFD 
codes  such  as  STAR-CD.  The  ADH  model  was  further  tested  in  this  particular 
study  of  modeling  the  flows  in  Ice  Harbor  forebay  using  a  flow  condition  chosen 
to  simulate  one  that  has  been  successfully  modeled  by  the  IIHR.  This  flow  con¬ 
dition  was  run  so  that  results  with  the  IIHR  model  could  be  compared  to  the 
ADH  simulations.  The  pool  elevation,  tainter  gate  openings,  and  powerhouse 
loadings  drive  the  discharges  listed  in  Table  2.  In  the  IIHR  study  for  the  Walla 
Walla  District,  this  flow  condition  was  referred  to  as  the  Base  Condition  1. 


Table  2 

Flow  Conditions  for  Base  Condition  1,  Total  River  Discharge 
74,700  cfs,  Pool  el  439.5 


Powerhouse  Unit  (Left  to  Right 
Looking  D/S) 

Discharge, 

cfs 

Spillway  Bay  (Left  to  Right 
Looking  D/S) 

Discharge, 

cfs 

1 

10,000 

1 

4,300 

2 

0 

2 

5,150 

3 

10,000 

3 

4,300 

4 

0 

4 

4,300 

5 

0 

5 

4,300 

6 

10,000 

6 

4,300 

7 

4,300 

8 

4,300 

9 

5,150 

10 

4,300 

Boundary  Conditions 

The  boundary  conditions  employed  to  model  this  particular  flow  condition 
were  set  similarly  to  the  validation  conditions  described  previously  and  consisted 
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of  setting  flow  rates  at  each  of  the  powerhouse  intakes,  each  spillway  gate,  at  the 
river  cross  section  making  up  the  upstream  limit  of  the  model,  at  the  water  sur¬ 
face,  and  at  all  solid  boundaries  comprising  the  computational  domain.  Each 
spillway  bay  also  has  a  geometric  configuration  associated  with  the  opening  of 
the  bay’s  tainter  gate.  Boundary  conditions  at  the  solid  boundaries  prohibit  the 
flux  of  mass  and  momentum  through  the  boundary.  This  no-flux  boundary  con¬ 
dition  is  applied  to  the  solid  boundaries  and  a  drag  coefficient  is  assigned  to  the 
face  of  each  element  forming  the  solid  boundaries.  The  water  surface  is  modeled 
as  a  rigid  lid,  which  is  reasonable  for  large,  deep  bodies  of  water  where  the 
Froude  number  is  small,  such  as  the  case  of  the  flow  approaching  Ice  Harbor 
Dam.  A  stress  free  condition  is  applied  for  the  tangential  component  of  shear  at 
the  water  surface. 


Inflow  boundary  conditions 

Flux  boundaries  are  prescribed  as  shown  in  Table  2.  The  inflow  boundary  is 
the  sum  of  the  discharges  through  the  powerhouse  and  the  spillway.  Velocity 
distribution  at  the  inflow  boundary  was  not  uniform,  but  rather  specified  in  a 
distribution  similar  to  that  used  by  IIHR  in  the  U2RANS  model.  This  distribution 
was  based  on  field  data  provided  by  the  Walla  Walla  District. 


Outflow  boundary  conditions 

Outflow  discharges  from  the  hydraulic  features  listed  in  Table  2  were  speci¬ 
fied  as  natural  boundary  conditions  to  each  element  face  making  up  each  outflow 
boundary.  The  average  velocity  at  each  feature  was  applied  to  each  outflow  face. 


Eddy  Viscosity 

Description 

As  was  mentioned  in  the  Governing  Equations  section  of  this  report,  the 
RANS  equations  are  closed  using  the  Boussinesq  method  of  adding  an  eddy  vis¬ 
cosity  to  the  molecular  eddy  viscosity  in  the  momentum  equations.  This  addi¬ 
tional  term,  often  referred  to  as  the  Reynolds  stress,  is  a  means  of  accounting  for 
the  effects  of  turbulence.  Here,  the  Reynolds  stress  is 

du. 

(3) 


Here,  r#  is  stress  acting  in  the  direction  xt  on  a  face  perpendicular  to  Xj.  The 
kinematic  eddy  viscosity,  Vy,  can  be  written  in  tensor  form  for  a  Cartesian  coor¬ 
dinate  system  as 
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V  V 

xy  xi 


V  V  V 

yx  yy  yz 


V  V 

zy  zi 


The  tensor  is  symmetric  so  that 


Vxy  tyx 


Vxz=  Vzx 


Vyz=  Vzy- 

These  turbulent  viscosities  can  be  grouped  in  terms  of  horizontal  and  vertical 
momentum  diffusion  due  to  turbulence.  The  terms  vxx ,  and  v^can  be  consid¬ 

ered  the  horizontal  eddy  viscosity  terms  and  the  terms  vxz ,  vyz,  and  vzz  represent 
the  vertical  eddy  viscosity.  The  Ice  Harbor  3D  model  was  tested  as  to  its  sensi¬ 
tivity  to  the  choice  of  horizontal  and  vertical  viscosity  values. 


Viscosity  sensitivity 

The  magnitudes  of  eddy  viscosity  were  varied  in  a  systematic  process  to 
determine  the  model’s  sensitivity  to  selection  of  both  horizontal  and  vertical  eddy 
viscosity.  This  is  particularly  important  in  this  study  where  the  length  scales  in 
the  horizontal  and  vertical  are  so  different.  Although  the  vertical  scale  (flow 
depth)  is  large  for  typical  river  models  (Ice  Harbor  forebay  is  as  much  as  140  ft 
deep),  the  horizontal  length  is  much  larger.  The  diffusion  of  momentum  due  to 
turbulence  becomes  important  as  the  flow  approaches  the  dam,  especially  the 
accelerating  flow  into  the  intakes  of  the  powerhouse. 

Flow  solutions  were  obtained  using  four  combinations  of  values  for  the  hori¬ 
zontal  and  vertical  eddy  viscosities.  The  selection  of  viscosities  tested  is  provided 
in  Table  3. 


Table  3 

Eddy  Viscosity  Values  Tested,  Total  River  Discharge  74,700  cfs, 

Pool  el  439.5 

Case  Number 

Horizontal  Eddy  Viscosity,  ft2/sec 

Vertical  Eddy  Viscosity,  ft2/sec 

1 

5.0 

5.0 

2 

5.0 

0.5 

3 

5.0 

0.25 

4 

0.5 

0.25 
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Figures  53-55  show  velocity  magnitude  contours  on  horizontal  planes  at  10-, 
40-,  and  80-ft  depths.  Figures  56  and  57  show  contours  of  velocity  magnitude  on 
vertical  planes  passing  through  powerhouse  unit  6  and  spillway  bay  5,  respec¬ 
tively.  Figures  58  and  59  show  the  magnitude  of  acceleration  at  powerhouse  unit 
6  and  spillway  bay  5,  and  Figures  60  and  61  show  hydraulic  strain  rate  compari¬ 
sons  at  the  same  locations. 

Model  results  with  Case  1  conditions,  horizontal  and  vertical  eddy  viscosity 
values  set  equal  to  the  vertical  values  of  5.0  ft2/sec,  are  seen  on  the  contour  plots 
as  being  somewhat  “stiff’  in  that  there  is  inadequate  diffusion  of  the  vertical 
momentum  due  to  turbulence. 

The  next  run,  Case  2,  was  conducted  using  a  horizontal  value  of  5.0  ft2/sec, 
but  with  the  vertical  eddy  viscosity  reduced  an  order  of  magnitude  to  0.5  ft2/sec. 
These  simulation  results,  shown  in  the  contour  plots,  increased  the  vertical  varia¬ 
tion  in  velocity  as  can  be  seen  by  comparing  the  velocity  magnitudes  with  those 
for  the  larger  vertical  viscosity. 

The  Case  3  conditions  again  used  a  horizontal  viscosity  of  5.0  ft2/sec,  but 
reduced  the  vertical  viscosity  in  half  to  a  value  of  0.25  ft2/sec.  This  reduction  in 
vertical  viscosity  had  no  noticeable  effect  on  the  resulting  solution.  Comparison 
of  the  different  cases  shown  in  Figures  56  and  57  show  that  the  vertical  variations 
of  velocity  magnitudes  are  essentially  the  same  as  a  result  of  using  vertical  vis¬ 
cosity  of  0.5  ft2/sec  and  0.25  ft2/sec. 

The  sensitivity  of  the  computed  solution  to  variations  in  the  horizontal  vis¬ 
cosity  was  also  tested.  The  horizontal  viscosity  value  for  the  Case  4  conditions 
was  reduced  an  order  of  magnitude  to  0.5  ft2/sec  and  the  vertical  viscosity  was  set 
to  0.25  ft2/sec.  Figures  53-55  show  that  this  reduction  of  horizontal  viscosity  did 
not  significantly  influence  the  velocity  distribution  throughout  the  flow  domain. 
The  horizontal  eddy  viscosity  of  5.0  ft2/sec  produced  the  appropriate  diffusion  of 
the  horizontal  momentum  due  to  turbulence. 


Comparison  of  Results 

The  horizontal  eddy  viscosity  value  of  5.0  ft2/sec  and  a  vertical  eddy  viscos¬ 
ity  value  of  0.25  ft2/sec  are  reasonable  choices  for  modeling  forebay  flows  at 
projects  such  as  the  Ice  Harbor  Dam.  The  flow  domain  was  not  identical  for  the 
ADH  and  U2RANS  models.  The  difference  is  especially  pronounced  at  the  right 
bank  line  (looking  downstream).  The  ADH  mesh  closely  models  the  actual  water 
line.  This  is  apparent  when  the  right  descending  bank  line  shown  in  the  prototype 
photograph  on  Figure  3  is  compared  with  the  CAD  model  shown  in  Figure  4. 
Also,  details  around  the  structure  were  different  in  the  two  models.  The  U2RANS 
reproduced  the  upper  end  of  the  lock  structure  in  more  detail  than  the  ADH 
model.  Also,  the  U2RANS  model  had  fine  resolution  near  the  structure’s  surface. 
This  will  be  evident  later  in  the  comparison  of  the  models’  velocity  gradients. 

Secondary  flow  variables  were  also  computed.  Derived  variables,  such  as 
acceleration  and  strain  rate,  are  used  by  other  researchers  in  modeling  migratory 
fish  behavior.  So,  these  flow  variables  can  be  important  if  the  computed  flow 
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field  is  used  to  forecast  fish  passage.  The  flow  acceleration  is  a  vector  quantity 
that  at  steady  state  can  be  expressed  as: 

V-VV  (5) 

The  hydraulic  strain  rate  is  here  defined  as  a  scalar  magnitude  of  the  strain- 
rate  tensor 


8ui 


dx. 


or  in  a  Cartesian  coordinate  system  as 


(6) 


du 

du 

du 

dx 

¥ 

dz 

8v 

dv 

dv 

dx 

dy 

dz 

dw 

dw 

dw 

dx 

dy 

The  comparison  of  ADH  and  U2RANS  results  are  shown  in  Figures  62-74. 
Figures  62-67  are  contours  of  velocity  magnitudes  at  10-,  40-,  and  80-ft  depths. 
Figures  68  and  69  are  velocity  magnitude  contours  on  vertical  slices  (plane  nor¬ 
mal  to  dam  axis)  at  powerhouse  unit  6  and  spillway  bay  5,  respectively. 

Figures  70-72  show  strain  rate  contours  at  10-,  40-,  and  80-ft  depth.  Figures  73 
and  74  show  strain  rate  contours  at  powerhouse  unit  6  and  spillway  bay  5.  The 
most  noticeable  differences  between  the  results  are  seen  in  the  section  views 
normal  to  the  spillway  and  powerhouse  although  these  differences  in  magnitudes 
are  small.  The  solutions  are  similar  and  the  use  of  either  to  make  engineering 
decisions  would  result  in  the  same  conclusions. 
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5  Summary  and  Conclusion 


Two  test  cases  were  made  a  part  of  this  ADH  model  validation  study.  The 
first  case  was  the  validation  of  the  ADH  model  by  comparison  to  observed  field 
data  at  the  Ice  Harbor  Dam  forebay.  This  tested  the  model’s  ability  to  simulate 
forebay  conditions  given  project  operations.  Comparison  of  the  model  and  field 
data  showed  the  numerical  model  was  capable  of  reproducing  field  conditions. 
This  first  test  also  served  to  evaluate  mesh  resolution  issues. 

The  second  test  was  directed  at  answering  the  question  of  how  well  the 
unstructured  ADH  solutions  matched  the  solutions  previously  computed  using 
the  IIHR-Hydroscience  and  Engineering’s  U2RANS  model.  The  particular  issue 
to  be  addressed  was  in  regards  to  the  computation  of  velocity  gradients  as  well  as 
the  primary  flow  variables  of  velocity  and  pressure.  Extensive  testing  to  ADH 
results’  sensitivity  to  choice  of  eddy  viscosity  was  also  included  in  this  part  of 
the  study.  The  solutions  from  each  model  were  very  similar. 

This  study  has  demonstrated  the  capability  of  modeling  the  approach  flow  to 
dams.  The  Navier-Stokes  module  of  the  ERDC’s  ADH  model  is  a  viable  means 
of  simulating  the  flow  approaching  and  interacting  with  hydraulic  structures. 
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Figure  1 .  Vicinity  map  of  Ice  Harbor  Dam 


Figure  2.  Ice  Harbor  Lock  and  Dam  (forebay  at  top  of  photo) 


Figure  3.  Ice  Harbor  Lock  and  Dam  (powerhouse,  spillway,  and  navigation  lock) 


Figure  4.  Overall  view  of  CAD  model  of  Ice  Harbor  forebay  (water  surface 
removed) 


Figure  5.  CAD  model  of  Ice  Harbor  forebay 


Figure  6.  Bathymetric  features  of  Ice  Harbor  forebay 


Figure  7.  CAD  model  of  powerhouse  and  adjacent  spillway 


Figure  8.  Detail  view  of  CAD  model  of  powerhouse 


Figure  9.  Details  of  powerhouse  and  adjacent  spillway 


Figure  10.  CAD  model  of  right  end  of  spillway  and  lock  floating  guard  wall 


Figure  1 1 .  Overall  view  of  surface  mesh  of  water  surface 


Figure  12.  Surface  mesh  of  typical  spillway  bay  with  tainter  gates  and  piers 


Figure  13.  Detail  view  of  surface  mesh  of  powerhouse  units 


Figure  14.  Surface  mesh  of  powerhouse  and  adjacent  bathymetry 


Figure  15.  Surface  mesh  of  spillway  near  navigation  lock  and  floating  guard  wall 


Figure  16.  Surface  mesh  (triangles)  and  portion  of  volume  mesh  (tetrahedra) 
within  powerhouse  unit  6  intake 


Figure  17.  Surface  mesh  (triangles)  and  portion  of  volume  mesh  (tetrahedra)  at 
spillway  bay  5 
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Figure  18.  Velocity  magnitude  contours  at  10-ft  depth 
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Figure  19.  Velocity  magnitude  contours  near  the  dam  at  10-ft  depth 


Figure  20.  Velocity  magnitude  contours  approaching  the  powerhouse  at  10-ft  depth 


Figure  21 .  Velocity  magnitude  contours  approaching  the  spillway  at  10-ft  depth 
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Figure  22.  Velocity  magnitude  contours  at  20-ft  depth 
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Figure  23.  Velocity  magnitude  contours  near  the  dam  at  20-ft  depth 


Figure  24.  Velocity  magnitude  contours  approaching  the  powerhouse  at  20-ft  depth 


Figure  25.  Velocity  magnitude  contours  approaching  the  spillway  at  20-ft  depth 
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Figure  26.  Velocity  magnitude  contours  at  40-ft  depth 


Figure  27.  Velocity  magnitude  contours  near  the  dam  at  40-ft  depth 


Figure  28.  Velocity  magnitude  contours  approaching  the  powerhouse  at  40-ft  depth 


Figure  29.  Velocity  magnitude  contours  approaching  the  spillway  at  40-ft  depth 


Figure  30.  Velocity  magnitude  contours  at  60-ft  depth 


Figure  31 .  Velocity  magnitude  contours  near  the  dam  at  60-ft  depth 


Figure  32.  Velocity  magnitude  contours  approaching  the  powerhouse  at  60-ft  depth 


Figure  33.  Velocity  magnitude  contours  approaching  the  spillway  at  60-ft  depth 


Figure  34.  Velocity  magnitude  contours  at  80-ft  depth 


Figure  35.  Velocity  magnitude  contours  near  the  dam  at  80-ft  depth 


Figure  36.  Velocity  magnitude  contours  approaching  the  powerhouse  at  80-ft  depth 


Velocity  Mag  (ft/sec) 


I  1.4 

I  12 

I  1 


0.8 

0.6 

0.4 

■  0.2 
0 


Q(kcfs) 

River  Total: 

105.5 

Powerhouse(1,2,4,6):  65.5 

Spillway: 

W.S.  439.5 

40.0 

ih_vaNd_Vmag_80-ft_swlay 


Figure  37.  Velocity  magnitude  contours  approaching  the  spillway  at  80-ft  depth 
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Figure  38.  Velocity  magnitude  contours  (on  a  vertical  plane  normal  to  dam  axis)  of  flow  entering 
center  intake  of  powerhouse  unit  6 
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Figure  39.  Velocity  magnitude  contours  (on  a  vertical  plane  normal  to  dam  axis)  of  flow  through 
the  center  of  spillway  bay  5 


Figure  40.  Flow  paths  colored  according  to  velocity  magnitude 


Figure  41 .  Close  view  of  flow  distribution  across  the  dam 
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Figure  42.  Comparison  of  prototype  data  (red)  and  CFD  results  (black)  at  10-ft  depth 
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Figure  43.  Comparison  of  prototype  data  (red)  and  CFD  results  (black)  at  20-ft  depth 
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Figure  44.  Comparison  of  prototype  data  (red)  and  CFD  results  (black)  at  40-ft  depth 
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Figure  45.  Comparison  of  prototype  data  (red)  and  CFD  results  (black)  at  60-ft  depth 
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Figure  46.  Comparison  of  prototype  data  (red)  and  CFD  results  (black)  at  80-ft  depth 
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Figure  47.  Prototype  data  at  20-ft 
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Figure  48.  Comparison  of  unrefined  mesh  (153,602  nodes)  and  refined  mesh  (548,464  nodes)  CFD 
results  at  10-ft  depth 


Figure  49.  Comparison  of  unrefined  mesh  (153,602  nodes)  and  refined  mesh  (548,464  nodes)  CFD 
results  at  40-ft  depth 
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Figure  50.  Comparison  of  unrefined  mesh  (153,602  nodes)  and  refined  mesh  (548,464  nodes)  CFD 
results  at  80-ft  depth 
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Figure  51 .  Comparison  of  unrefined  mesh  (153,602  nodes)  and  refined  mesh  (548,464  nodes)  CFD 
results  at  powerhouse  unit  6 
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Figure  52.  Comparison  of  unrefined  mesh  (153,602  nodes)  and  refined  mesh  (548,464  nodes)  CFD 
results  at  spillway  bay  5 
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Figure  53.  Velocity  magnitude  contours  of  flow  at  1 0-ft  depth 
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Figure  54.  Velocity  magnitude  contours  of  flow  at  40-ft  depth 
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Figure  55.  Velocity  magnitude  contours  of  flow  at  80-ft  depth 
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Figure  56.  Velocity  magnitude  contours  of  flow  entering  powerhouse  unit  6 
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Figure  57.  Velocity  magnitude  contours  of  flow  through  spillway  bay  5 
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Figure  58.  Acceleration  magnitude  contours  of  flow  entering  powerhouse  unit  6 
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Figure  59.  Acceleration  magnitude  contours  of  flow  through  spillway  bay  5 
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Figure  60.  Strain  rate  contours  of  flow  entering  powerhouse  unit  6 
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Figure  61 .  Strain  rate  contours  of  flow  through  spillway  bay  5 
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Figure  62.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  at  10-ft 
depth 
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Figure  63.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  near  the 
spillway  at  10-ft  depth 
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Figure  64.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  at  40-ft 
depth 
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Figure  65.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  near  the 
spillway  at  40-ft  depth 
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Figure  66.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  at  80-ft 
depth 


0  100  200  300 

U2RANSBC1 


8100 


8000 


7900 


7800 


7700 


£ 

>\—i 

200 

X.ft 


300  400  500 

ADH  BC1 


BC1 

River  74.7kCfS,  PM:  30kCfS,  SWA4.7kCfS,  WSEL  439.5 
Ph :  units  1,3,6:  lOkcrseach 

SW:  Bays  1,3,4,5,67,8,1 0: 4.3kcfs  each,  Bays2,9: 5.15kcfs  each 


I 


Velocity  Mag  (ft/sec):  0  0.2  0  4  0.6  0.8  1  1.2  1.4  1.6  1.8  2 


ih_cfti5005_vs_u2rans_80-lt_ph  lay 


Figure  67.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  near  the 
powerhouse  at  80-ft  depth 
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Figure  68.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  (on  a 
vertical  plane  normal  to  dam  axis)  of  flow  entering  powerhouse  unit  6 


Figure  69.  U2RANS  (left)  and  ADH  (right)  model  results  showing  velocity  magnitude  contours  (on  a 
vertical  plane  normal  to  dam  axis)  of  flow  entering  spillway  bay  5 
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Figure  70.  U2RANS  (left)  and  ADH  (right)  model  results  showing  strain  rate  contours  near  the  spillway  at 
1 0-ft  depth 
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Figure  71 .  U2RANS  (left)  and  ADH  (right)  model  results  showing  strain  rate  contours  in  the  forebay  at 
40-ft  depth 
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Figure  72.  U2RANS  (left)  and  ADH  (right)  model  results  showing  strain  rate  contours  near  the 
powerhouse  at  80-ft  depth 
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Figure  73.  U2RANS  (left)  and  ADH  (right)  model  results  showing  strain  rate  contours  (on  a  vertical  plane 
normal  to  dam  axis)  of  flow  entering  powerhouse  unit  6 
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Figure  74.  U2RANS  (left)  and  ADH  (right)  model  results  showing  strain  rate  contours  (on  a  vertical  plane 
normal  to  dam  axis)  of  flow  entering  spillway  bay  5 


